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Equation of state of a meson gas 
from the PNJL model for interacting quarks 

Kanako Yamazaki* and T. Matsur 
Institute of Physics, University of Tokyo, Komaba, Tokyo, Japan 

We compute the equation of state of hadronic matter at finite temperature with zero net baryon 
density by the Nambu-Jona-Lasinio model for interacting quarks in uniform background temporal 
color gauge fields. In the low temperature confining phase, unphysical thermal quark-antiquark 
excitations which would appear in the mean field approximation, are eliminated by enforcing van- 
ishing expectation value of the " Polyakov-loop" of the background gauge field, while in the high 
temperature confining phase this expectation value is taken as close to unity allowing thermal ex- 
citations of free quarks and antiquarks. Mesonic excitations in the low temperature phase appear 
in the correlation energy as contributions of collective excitations. We describe them in terms of 
thermal fluctuations of auxiliary fields in one-loop (Gaussian) approximation, where pions appear 
as Nambu-Goldstone modes associated with dynamical symmetry breaking of the chiral symmetry 
in the limit of vanishing bare quark masses. We show that the equations of state reduces to that of 
free meson gas with numerically small corrections arising from the composite nature of mesons. 

PACS numbers: 11.10. Wx,11.30.Rd, 12.38. Mh,25. 75. Nq 



The Nambu-Jona-Lasinio (NJL) model, originally in- 
troduced to describe profound consequences of presumed 
underlying chiral symmetry in hadron dynamics in anal- 
ogy to the BCS theory of superconductivity has been 
widely used to describe effective quark dynamics respect- 
ing the symmetry . When it is applied in the mean field 
approximation to a finite temperature system, however, 
it generates thermal excitations of quark quasiparticles 
even at low temperatures, where quarks are to be con- 
fined in individual mesons and baryons. 

A phenomenological model which removes these 
unphysical modes of excitations was proposed by 
FukushimaQ by making use of the Polyakov loop which 
is known to characterize confinement-deconfinement 
transition in strong coupling pure gauge theories. [4j This 
procedure eliminate the quark single particle distribu- 
tion by phase cancellations, leaving behind excitations 
of triad of quark quasiparticles in color singlet con- 
figuration. The quark-liberating deconfining transition 
was described as a change of the expectation value of 
the Polyakov loop treated as an order parameter with 
effective Landau- Ginzburg type free energy functional. 
Fukushima's model has been reformulated as a mean field 
theory in a uniform background temporal gauge field 
and has been studied extensively by other groups |5|- 
9]. Their model, referred to as the PNJL model, allows 
to incorporate mesonic correlations beyond mean field 
approximation Udcli as had been studied previously in 
the original NJL modelfllL . in analogy to the classic 
example of plasma oscillation. (l3j Numerical computa- 
tion with the model has been shown to reproduce the 



equation of state calculated by the Monte Carlo simula- 
tion of the lattice gauge theory at small non-zero chemical 
potential rather well. The model has been used to com- 
pute susceptibilities of hadronic matter at finite chemical 
potential 14] , studied in the regime of complex chemical 
potential 15[ and very recently used to study the effect 
of confinement in the study of QCD phase diagramflfjj. 

It is the purpose of these notes to re-examine the model 
at zero chemical potential and show that equation of 
state of a meson gas can be derived explicitly at low tem- 
peratures by the method of auxiliary fields which physi- 
cally express effective meson fields build-up as a quark- 
antiquark bound states. Similar results have been ob- 
tained earlier by others by numerical computations of 
the mesonic correlation energies with slightly different 
models 0, H|, but we will show explicitly that in the chi- 
ral limit the equation of state becomes that of a gas of 
mesons and a correction due to the composite nature of 
mesons as the individual excitations of the quark triad. 
It is found that the former dominates the equation of 
state at low temperatures and the latter do not alter 
the free meson gas result significantly. As the temper- 
ature increases the mean field contribution from under- 
lying quark quasiparticle excitations become important 
due to the non- vanishing Polyakov loop expectation value 
and will eventually dominate the free energy in deconfin- 
ing phase. Collective mesonic excitations melt into the 
continuum of quark anti-quark excitations. 

We begin with an Euclidean path-integral expression 
of the partition function for the NJL model of quarks in 
external temporal color gauge field: 



Z(T,A 4 ) = J[dq][dq]exp J d 4 a;>C N j L (g, q, A4) 



(1) 
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where we have used an abbreviate notation for integral 
over the space and imaginary time (drd 3 r — d 4 x), from 
— oo to oo for each of three space coordinates and from 



2 



to f3 = l/T for imaginary time r, and the effective 
Lagrangian is given by 

£njl(<7, q,M) = q{i^D p ~m )q + G [(qq) 2 + {iqj 5 rq) 2 ] 

(2) 

for two flavor light quarks, q = (it, d) in external (clas- 
sical) temporal color gauge field, = + gAoS^fi 
and mo is the bare quark mass which breaks the chi- 
ral symmetry explicitly. Here we have used a standard 
Minkovsky metric notation with the real time being re- 
placed by the imaginary time r = ixo and the time com- 
ponent of the gauge field replaced by A4 = iAq [l7l |. 
We note here that quark fields are the only dynamical 
variables which describes thermal excitations of the sys- 
tem, while the temporal component of SU(3) gauge fields 
A 4 = \\ a Al, \ a being the 3 x 3 Gell-Mann matrices, 
only plays a side role imposing constraints on color con- 
figurations of thermal quark excitations. In the following 
calculation, we take Ansatz of diagonal representations 
for A4 along with the earlier works [3, [B-Q ■ 

To describe mesonic excitations, we introduce four 
auxiliary bosonic fields 4>i — (a, tt) coupled to quark 
densities (qq,iq^/5Tq) by multiplying Z(T, A4) by a con- 
stant Gaussian integral. This procedure, known as the 
Hubbard-Stratonovich transformation, converts the four 
point NJL quark interaction to Yukawa coupling of the 
bosonic fields 4>i to corresponding four quark densities. 
Performing the Grassmann integral over the quark Dirac 
fields, the partition function becomes 



Z{T,A 4 ) = J[d<t>]< 



-i(4>Ai) 



(3) 



where the exponent is given by the integral 

i"0M4) = ^ J d A x ((a - m ) 2 + tt 2 ) 

-Trln [j3 (i^D M + a + i Tj t ■ tt)} (4) 

where the trace in the second term is taken over the ar- 
guments of quark fields, including the space-time coordi- 
nates with anti-periodic boundary condition in the imag- 
inary time axis, Dirac gamma matrices and the isospin 
and color indices. 

We evaluate the functional integral over the mesonic 
auxiliary fields by the method of steepest descent. Let 
<po(ao,7To) give local minimum value of integrand. Here- 
after we choose itq = 0, along with 00 = Mo, which 
evidently satisfy stationary conditions for effective pion 
fields. Expanding the effective action in a power series 
of <j) — 4>q up to the quadratic term and performing the 
gaussian integral over the bosonic fields in eq. © , we find 



Q(T,A 4 )=T U 



iTr M ln- 
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(5) 



up to one-loop fluctuations of the collective meson fields, 
where the trace is to be performed over the internal co- 
ordinates of the auxiliary meson fields. 



The thermodynamic potential in the mean field ap- 
proximation, £Imf(T, A4), or the corresponding pressure 
Pmf(T, A4), is retrieved from the leading term Iq by the 
relation: D,mf(T, A4) = TIq — —pmfV. Converting the 
coordinate representation to the momentum representa- 
tion and evaluating the discrete sum over the quark Mat- 
subara frequencies by the standard method of contour 
integration|18j. we find 



p MF (T, M) = 2N f J ^3 ^tr c [f(E p + igJU) 

+ f(E p -igA 4 )}+p MF (M Q ) (6) 

where the first term is the pressure of the quark quasi- 
particle in the external gauge fields with distribution, 



f(E p ± igM) 



1 



e li(E p ±igAi) _|_ 1 



(T) 



with energy E p = ^Jp 2 + Mq and 

' A d 3 p 



f d A v I 

P° MF (Mo) = 3 x 2N f J j^E p - —(Mo - m ) 2 

(8) 

is the pressure exerted by the quark quasiparticles in the 
"Dirac sea": It requires momentum cut-off A. The fac- 
tor 3 in the first term of (JSJ accounts for color and the 
factor 2 for spin degeneracies. Note that this "vacuum 
pressure" does not depend on the external background 
gauge potential due to the cancellation of the effects of 
the potential on particle and antiparticle. It depends 
on temperature indirectly through the temperature de- 
pendence of the effective quark mass Mo. The vacuum 
pressure should be removed by the renormalization con- 
dition. The quark pressure should also be removed in 
the low temperature phase where quarks are confined in 
hadrons. 

In the mean field approximation, the value of the quark 
quasiparticle mass M is determined by the stationary 
condition of Iq with respect to Mo. This procedure leads 
to 



M - m = 8GN f tT c 



d 3 p M 



(2irf E p 
- f(E p - igJU)) 



{I- f{E p + igA i ) 
(9) 



In the above expressions, the constant temporal gauge 
field Ai appears as a phase factor together with the quark 
quasiparticle energy in the single particle distribution 
function ([7]). It looks very similar to the gauge invari- 
ant (path-ordered) Polyakov loop phase integral, 



L(r) = "Pexp 



dr Ai(r, r) 



(10) 



whose thermal expectation value measures the extra free 
energy associated with the color charge in fundamental 
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representation fixed at a spatial point r. Although this 
connection is not strict, since quarks are moving in a 
uniform background gauge field not fluctuating either in 
space or in imaginary time, we adopt the procedure of 
Fukushima, along with the TUM group, to replace the 
phase factor in the quark quasiparticle distribution func- 
tion by the thermal average of the Polyakov loop. 
We replace {\tv c f{E p + igA A )) by [ijj 

, . <le 2 ^+2$e^ + l 

U{E p ) - e3fjEp + 3 ^ e 20E p + 3$ e 0B p + 1 ^ > 

where $ = -|(tr c L), and $ = ■|(tr c L'!'}. The anti-quark 
distribution, f<j>(E p ) = (^tr c f(E p — igA^)), is obtained 
from the quark distribution by the replacement of $ by 
$ and vice versa. At zero chemical potential, two distri- 
butions are identical and <£■ is real so that $ = $. Very 
interesting observation here is that although in the de- 
confining phase where $ = $ = 1 , one would simply have 
ordinary quark distribution, 

/*(£p)U=i = &e ,E v + l > ( 12 ) 

in the confining phase, where $ = $ = 0, we have in- 
stead, 

/*CEp)l*=0 = e 3(SE p + l ( 13 ) 

which may be interpreted as a distribution function of 
triad of quark quasiparticles exciting together in color 
singlet configuration. It looks something like a baryon 
although no effect of interaction is taken into account 
between three quarks to form a baryon, as reflected by 
the same value of momenta for all three quarks. We 
also note that the degrees of freedom of excitations is 
not reduced to 1/3 of color triplet free quark excitations. 
However, there still exist reduction of effective degrees of 
freedom due to the change of the phase space of quark 
momentum p and that of quark triad momentum ip: The 
"triad" carries the energy E tr i = 3E p , the momentum 
Ptri = 3p and the mass M tr i = 3 Mo. This implies that 
the Lorentz invariant phase space integral over the quark 
momentum is reduced when written in term of the triad 
momentum. 

Here we determine $ phenomenologically by adding 
an effective potential U{T,<&): fi(T, $) = (fl(T,A 4 )} + 
VU(T, $). The parameters of U(T, $) are chosen so that 
$ = for the low temperature confining regime where 
thermal quark quasi-particle excitations are forbidden, 
while $ ~ 1 above deconfining temperature, quite op- 
posite to the normal behavior of the usual order param- 
eter. We note here that, although gluon fields are not 
treated as dynamical variables, gluon excitations may be 
included in a phenomenological fashion through the ef- 
fective potential by setting, U(T, $ = 1) = -^T 4 so 
that the equation of state matches to that of a free quark- 
gluon plasma at asymptotically high temperatures. 



The equation of state in the mean-field approximation 
to the PNJL model is now given as a maximum of 

+Pu F (M ) - Ap vac (14) 

with respect to variation of $ and the quark quasipar- 
ticle mass Mo, which measures the magnitude of chiral 
condensate, is determined by 

M -mo = l2GN f J (l - h(E p ) - U(E p )) . 

(15) 

The constant Ap vac is chosen so that the pressure be- 
comes zero at zero temperature and the cut-off at A in 
the momentum integral in (|15l) is applied only to the 
vacuum component. 

We show the temperature dependence of <E> and M ob- 
tained from these mean- field equations in Fig. 1. In the 
chiral limit m = 0, the gap equation (IT5|) possesses a 
non-trivial solution Mo only at temperatures below 
a critical temperature T c . In the ordinary NJL model, 
with no "Polyakov loop prescription", the chiral tran- 
sition takes place at relatively low temperature, while 
in the pure gauge lattice simulation, the confinement- 
deconfinement transition takes place at higher temper- 
ature as a first order transition with the discontinuous 
change in the value of $. In the PNJL model, these 
two distinct transitions interfere each other through the 
quark quasiparticle loops. As a result, the two transitions 
takes place at a similar temperature: the transition tem- 
perature for deconfinement becomes lower and becomes 
smooth crossover transition, while the chiral transition 
takes place at higher temperature, since the quark exci- 
tations are suppressed at low temperatures by small ex- 
pectation value of the Polyakov loop. With a finite bare 
quark mass mo, the chiral symmetry is broken explicitly 
and the chiral transition also become a smooth cross over 
transition. 

Now we come to our main task to evaluate mesonic 
correlation energy. The contribution of the mesonic cor- 
relation energy to free energy, J7m, or the corresponding 
pressure pm = —£Im/V> is given by 

Pu{T,A i ) = -^ J (0f n ^ ((2G9" 1 - n CT )] 

+31n[/3 2 ((2G)- 1 -n,)]} (16) 

where the sum is over the bosonic Matsubara frequencies 
Lj n = 2nirT which appear in the arguments of meson self 
energies: 

u a/ ^{w n , q ,A A ) = nl /7T (A A ) + nl /7T (u n ,q,A 4 )(17) 
The "tad pole" term 

n 1 -H 1 --2TVtr f-^P- - 

(18) 
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FIG. 1: Temperature dependence of effective quark mass Afo 
(red) and expectation value of the Polyakov loop (blue) : the 
dashed curves are for vanishing bare quark mass (the chiral 
limit) and the solid curves are with finite bare quark mass mo 
which is chosen to reproduce the pion mass m n — 140MeV. 



does not depend on w„ but the dispersive term 

I^K.g,^) - K + q 2 +4Af 2 )FK,q,^ 4 )(19) 
Ul(Lj n ,q,A A ) = (io 2 n + q 2 )FK,q,A 4 ) (20) 

contains explicit uj n dependence and 

d 3 p 1 



F(oj n ,q, A 4 ) 



(2tt)3 ( €n ,-gA 4 y + E2 
1 



{e n , - gA A + ui n y + E 



(21) 



p+q 



where e n = (2n + l)nT are the quark Matsubara frequen- 
cies. 

We now show that the contributions of the mesonic 
correlation to the pressure (|T6|) indeed contain those of 
free meson gas composed of massless pions and massive 
sigma mesons. For this purpose we first eliminate the 
"tadpole" term in ±1* and II 4 by using the stationary 
condition, or the "gap equation" (|9]) which determines 
the quark quasiparticle mass Mq in the symmetry broken 
phase. We then find, 



(2G)- 



(w 2 + q 2 + 4Af 2 )F(w„, q, M) + 



(2G)" 1 -n 7r = (ujl + q 2 )F(LJ n ,q,A i ) + 



m 

2GMn 



m 
2GM 
(22) 

(23) 



In the chiral limit, m Q — 0, we can thus separate the 
contributions of the collective bare meson modes from 
non-collective individual excitations 



p M (T, A,) = p c ^(T) +pTr coU iT,A 4 ) 



-coll. 



(24) 



The first term, independent of the color gauge potential, 
are identical to the pressures of free sigma meson gas 



with mass m a = 2Mq and of free massless pion gas with 
3 isospin degeneracy. 



pT{t) 



d 3 q 



2 

6u>q 3 



(25) 



where uj p = ^/q 2 + 4 M5 2 is the sigma meson energy and 
/s(w) = l/(e /3w — 1) is a bosonic single particle dis- 
tribution function. The mesonic vacuum pressure p° M 
which depends on Mq will be removed at zero tempera- 
ture by the renormalization, but we retain the deviation 
due to the shift of Mq at finite temperature. The mo- 
mentum integral for the pressure of the massless pion 
gas can be evaluated analytically and one finds a famil- 
iar Stephan-Boltzmann pressure with three-fold isospin 
degeracy: p% c n = 3 x f^T 4 = 0.33 x T 4 . 

The non-collective contribution depends on the color 
fields since it is composed of the underlying quark degrees 
of freedom. The sum over the bosonic sum over w„ can 
be performed again by contour integration and we find, 



non— col. 

Pm 



= -2 



du 



d 3 q 1 
(2tt) 3 2ni 7 

F(u + ie, q, A A ) 
F(uj - ie,q,Ai) 



xln 



(26) 



where F(u>, q, A4) consists of two terms 



F(w,q,A4.) = F scat (uj,q,A 4 ) + F pair (u;, q, A 4 ) (27) 
The first term is the scattering term given by 



F a . 



1 



d 3 p 



E„ — E, 



'p+q 



2 J (2tt) 3 E p E p+q (E p - E p+q y - ^ 
xtr c (f(E p - igAi) - f{E p+q - igA 4 )) 



(28) 

while the second term corresponds to the pair creation 
and annihilation term 



F ■ — 

± pair — 



d 3 P 



1 



E v + E; 



p+q 



1 

2 J (2tt) 3 E p E p+q {E p + E p+q y 2 - , 

xtr c (1 - f(E p - igAi) - f(E p+q - igM)) 

(29) 

Again, the external gauge fields appear as a phase factor 
in the quark distribution function as in the mean-field 
approximation. We therefore replace these phase factors 
by the Polyakov loops and then substitute them by sta- 
tistical average as in (fTTj) : 

F{uJ,q) = {F{u,q,Ai)) = J 7 aoat(w,g') + -Fpairv^g) 

where 7 r sca t(a;, q) and J 7 pa ir(w,g) are given by ([28)) and 
(|29[) . respectively, with the quark distribution replaced 
by CLE]). 
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The function T(ui ± ie, q) contains a real part T\ (w, q) 
and an imaginary part iJ-^w,*?). It is convenient to 
write 



where the argument cf> is given by 



(w, q) = tan 



Fi{u,q) 



(30) 



(31) 



The real part and imaginary part of the function J- are 
further decomposed into two parts. The two components 
of the real part are given by the principal part of the 
integrals, while the two components of the imaginary part 
contains the energy conserving <5-functions: 



d 3 p 



1 



{U{E P ) - f$(E p+q )) 



pair, 2 



A J {2^fE p E p+q 

x (8{uj + E p - E p+q ) - S(uj - E p + E p+q )) 

\ J (tr EpEp+q (1 ~ MEp) -M E ^ 

x (5(u + E P + E p+q ) - 5(uj - E p - E p+q )) 



It is evident that the scattering term has a non-zero 
imaginary part in the space-like energy-momentum re- 
gion (q 2 > w 2 ), while the pair creation/annihilation term 
is non- vanishing only in the time-like region (q 2 < u> 2 ). 
It is important to note that non-collective mesonic corre- 
lation arises only from non-vanishing imaginary part of 
T(q,ui). The pressure arising from the non-collective or 
individual excitations of the system is given by 



Pm 



i— col. 



(T) = -4 



dh. 



duj 



1 



(27T)- JO 
+ Ap non-col. 



2tt - 1 



2cj>(Lu,q) 



(32) 



where the temperature independent background pressure 
^P™"™' contains the vacuum pressure due to the fluc- 
tuation of the pair of quark-triad and anti-quark triad 
which is removed again when we estimate the physical 
pressure at finite temperature. 

We note that foregoing analysis only applies to the chi- 
ral limit (mo = 0) and symmetry broken phase. For non- 
zero value of mo, the separation of the collective mode 
and the individual excitation is not as simple as above 
analysis due to the term mo/(2GAio) in the dispersions 
(|2"21 and (|2U)) of each kind of mesons. Also at high tem- 
peratures we have no broken-symmetry solution in the 
gap equation. In such cases we have to go back to the 
original expression (|16[) in order to compute the mesonic 
correlation energy as has been done in [8j . Collective me- 
son poles appear separated from the continuum of indi- 
vidual excitations of a pair of quark triad and anti-quark 
triad and the former contribution dominates the pressure 
at low temperatures, yielding meson gas equation of state 
again. More detail discussion of this analysis is reported 
in EH. 



We show in Fig. 2 the pressure calculated by the 
present method. The upper panel (a) is the result com- 
puted in the chiral limit mo = 0. At low temperatures 
below the second order chiral transition temperature T c , 
the pressure from quark quasiparticle excitations are sup- 
pressed by the quenching of the distribution function: 
only excitations of triad of quarks with color singlet con- 
figuration are allowed with effective excitation energy 
three times that of massive quark excitation, this contri- 
bution to the pressure is hence suppressed. The dominant 
mode of excitations which determine the pressure is due 
to massless pion excitation with 3 isospin degeneracy as 
shown in the red solid curve. We found the contribution 
from massive sigma meson excitations and non-collective 
individual excitations of the quark triads and antiquark 
triads are also negligible at low temperatures. At high 
temperatures, massless quark excitations in the mean- 
field approximation dominate the pressure with compa- 
rable contributions of the gluon pressure phenomenolog- 
ically introduced in the construction of the effective po- 
tential. Contribution from the mesonic correlations be- 
comes negligible at high temperatures, and the pressure 
approaches to that of ideal massless quarks and gluons in 
our model. The deviation is due to the slight deviation of 
$ from unity. With symmetry breaking finite bare quark 
mass, shown in the lower panel (b), the pion becomes 
massive and its pressure decreases exponentially at very 
low temperature, otherwise the behavior of pressure at 
higher temperatures is qualitatively unchanged. 

We have shown that an equation of state of meson gas 
can be derived from interacting quarks with a Nambu- 
Jona-Lasinio type interaction in uniform background 
color field. In the low temperature confining phase, the 
unwanted quark quasiparticle excitations in the mean- 
field approximation are removed by Fukushima's method, 
introducing the quenched quark distribution in the van- 
ishing Polyakov-loop expectation value, which converts 
the distribution of individual quarks into that of triads 
of quark in color singlet configuration. The mesonic ex- 
citations which appear in the correlation energy as a 
correction to the mean-field approximation, dominates 
in the low temperature confining phase. We also found 
that the corrections due to non-collective individual ex- 
citation of triads of quarks or antiquarks are suppressed 
compared to collective meson excitations. Although each 
of these quark triads somewhat resembles a baryon car- 
rying unit baryon number; these quarks carry the same 
amount of energy-momentum, hence not correlated spa- 
tially. It remains to be seen how these quark-triads are 
converted to real baryons taking into account internal 
correlations among three quarks. For extensions of the 
present method to systems at finite baryon chemical po- 
tentials, such more realistic treatment of baryonic corre- 
lations would therefore be very important. 
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FIG. 2: Pressure scaled by T 4 as a function of tempera- 
ture: the upper panel (a) is for vanishing bare quark mass 
(the chiral limit) and the lower panel (b) is with finite bare 
quark mass mo which is chosen to reproduce the pion mass 
= 140MeV. In comparison, the pressures calculated in 
the mean-field approximation are shown with (without) gluon 
contributions. 
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